RELATIVISTIC FLUIDS AND THE 
PHYSICS OF GRAVITATIONAL 
COLLAPSE 



Questions 

Why is the study of gravitational collapse important 
for General Relativity? 

Why is the gravitational collapse a highly dissipative 
process. 

Is there any life between quasi-equilibrium and non- 
equilibrium? 

What happens if we relax the Pascal principle? (and 
why?) 

How does the electric charge and dissipation affects 
the evolution of massive stars? 



2 Why is the gravitational collapse a highly 
dissipative process? 
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3 Non— comoving (Bondi 1964). 
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3.1 Equilibrium a; = 0, and quasi— equilibrium u 2 « 
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v « 0. 



<j « 0(\\u,vu). (20) 



0(cj 2 ) = A 2 = z> 2 = Az> = A = v = 0, (21) 

3.2 Effective variables and the post— quasi— static ap- 
proximation: The life between quasi— equilibrium 
and non— equilibrium 
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3.3 The algorithm for modelling spheres out of equi- 
librium 

1. Take an interior (seed) solution to Einstein equations, 
representing a fluid distribution of matter in equilib- 
rium, with a given 

p st = p(r) P rst = P r {r) 

2. Assume that the r-dependence of P and p is the 
same as that of P rst and p st , respectively. 

3. Using equations (25) and (24), with the r dependence 
of P and p, one gets m and v up to some functions 
of t, which will be specified below. 

4. For these functions of t one has three ordinary dif- 
ferential equations (hereafter referred to as surface 
equations), namely: 

(a) uj = rs e^~ u ^ 2 evaluated on r = r^. 

(b) Tftp = evaluated on r = r^. 

(c) The equation relating the total mass loss rate 
with the energy flux through the boundary sur- 
face. 

5. Depending on the kind of matter under considera- 
tion, the system of surface equations described above 
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may be closed with the additional information pro- 
vided by the transport equation and/or the equa- 
tion of state for the anisotropic pressure and/or ad- 
ditional information about some of the physical vari- 
ables evaluated on the boundary surface (e.g. the 
luminosity) . 

6. Once the system of surface equations is closed, it may 
be integrated for any particular initial data. 

7. Feeding back the result of integration in the expres- 
sions for m and v, these two functions are completely 
determined. 

8. With the input from the point 7 above, and using 
field equations, together with the equations of state 
and/or transport equation, all physical variables may 
be found for any piece of matter distribution. 
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4 Relaxing the Pascal principle? (and why?) 

m + AnP r r 3 , ^ x 2 (Pj_ - P r ) 
r (r — 2mJ r 

Pr+P± 
2 

• What might be the origin of local anisotropy? 

1. "Exotic" phase transition (e. g. pion condensate). 

2. Magnetic fields 

3. Type II superconductor 

4. Type P superfluid 

5. Boson stars 

6. Viscosity 

7. Anisotropic velocity distributions 

8. Two fluid systems 

• How does the properties of the locally anisotropic 
system differs from the locally isotropic one? 

1. Cracking induced by perturbations of local isotropy 
of pressure. ( p± ~ Pr ) << 1 

2. Changes in the total mass allowed for a compact ob- 
ject 



5 Comoving (Misner and Sharp 1964). 

ds 2 _ = -A 2 dt 2 + B 2 dr 2 + {Cr) 2 {d6 2 + sin 2 6># 2 ), (27) 
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The electromagnetic energy tensor and the Maxwell 
equations 



E a(3 = ^ [FaF^ - -F l6 F l5 g aP 



(37) 



Fa/3 — 4>/3,a — 4>ajh 

F^-p = 4ttJ°, 

4> a = <$>5° a , J Q ' = ?V Q , 



(38) 
(39) 

(40) 



s(r) = 4tt f <;B{Crfdr, 



(41) 



.1' /;' a 2, , . , 



(42) 
(43) 



10 



sAB 

{Crf 



44) 



5.2 The Einstein equations 
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The exterior spacetime and junction conditions 



1 2M(v) | Q'' 
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6 Dynamical equations 

The non trivial components of the Bianchi identities, 
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U = tDtC = DtR < (in the case of collapse). 

(58) 
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^ + P r + 2e - ^o) P> T /7 = 
- ^ + P r + 2e - ^(j) m + 4tt (p r + e - ^770) P 3 - ^ 
^ (Pr + e - ^w) + 2(P r - Pl + e - 2rja) 
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Force = Mass (density) x Acceleration (65) 
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The transport equation 
1 Maxwell-Fourier law and causality 

q(x,t) = -KVT(x,t) 



de = 7<iT and 
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q(x, t + r) = -kVT{x , t) (75) 
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r « 10 11 sec. (phonon — elect.) 

r ~ 10~ 13 sec. (phonon — phonon) 

t « 10" 3 sec. #ei7 ~ 1.2° AT 

r « 1CT 1 — lCT 4 sec. neutron — star — matter 

(82) 



83) 



K T 2 g , t 

Uj 1 



2t 



3L7 1 1 



r " r£ 



1 



(84) 



18 



^ + P r + 2e - ^ar?) (1 - a)P> T /7 
= (l-a)F grav + F hyd + D R T 

T 

" T2q D T (^) - D T e 



+ E 



2t " * UT 2 . 
/5Z7 3 1 



1 

P 



(85) 



P 



grav 



= -(/i + P r + 2e-^7 



x 



/ 4 \ s 2 

m + 4tt ( P r + e - -r?cr I P 3 - — 



R 2 



,(86) 



+ 2(P r - Pj_ + e - 277a) 



1 



P 4ttP 4 



(87) 



a = — I [i + P r + ze — -(777 



88) 



Mm 



r](M + [P r ]+2[ e ]-|[ W ; 



x 10 



-42 



(89) 



19 



[k], [T], [t], en erg. s 1 cm 1 K l 1 K 1 s g. cm 

1. Pre-supernovae 

(Martinez 96) ([«] « 10 37 ; [T] « 10 13 ; [r] « 
[fi] « 10 12 ) a « 1. 

2. Inflation /i + P = a = 1 
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• Once the transport equation has been taken into ac- 
count, then the inertial energy density and the "pas- 
sive gravitational mass density" , i.e the factor multi- 
plying DtU and the first factor at the right of (64) re- 
spectively (which of course are the same, as expected 
from the equivalence principle), appear diminished 
by the factor 1 — a, a result already obtained, but 
here generalized by the inclusion of the viscosity and 
radiative phenomena. 

• As far as the right hand side of (85) is negative, the 
system keeps collapsing. However, let us assume that 
the collapsing sphere evolves in such a way that, for 
some region of the sphere, the value of a increases 
and approaches the critical value of 1. Then, as this 
process goes on, the ensuing decreasing of the gravi- 
tational force term would eventually lead to a change 
of the sign of the right hand side of (85). Since that 
would happen for small values of the effective inertial 
mass density, that would imply a strong bouncing of 
that part of the sphere, even for a small absolute 
value of the right hand side of (85). 

• Observe that the charge does not enter into the defi- 
nition of a. However it does affect the "active gravi- 
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tational mass" (the factor within the square bracket 
in (86)). 

• The repulsive Coulomb term (the last term in (87)) 
depends on Drs and always opposes gravitation. Its 
effect is reinforced if Drs is large enough to violate 
(67) , in which case the charge will decrease the "ac- 
tive gravitational mass" term in (86). 
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Conclusions 



Dissipative phenomena may play a relevant role in 
the dynamics of collapse. In particular, relaxational 
effects may drastically change the outcome of gravi- 
tational collapse. 

The dynamical regime may be approached by means 
of succesive approximations. 

Local anisotropy of pressure has to be taken into con- 
sideration in the study of the structure and evolution 
of massive stars. 

Not only electric charge but also its distribution may 
be relevant in stellar structure and evolution. 

Local (non-gravitational) effects may be crucial to 
determine the outcome of gravitational collapse. 

Learn to deal with comoving and non-comoving frames! 
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